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Abstract. The parallel code NMAGIC is an implementation of a particle-based 
method to create made-to-measure models in agreement with observations of galaxies. 
It works by slowly correcting the particle weights of an evolving N-body system, until 
a satisfactory compromise is achieved between the goodness of the fit to a given set of 
observational data, and some degree of smoothness (regularization) of the underlying 
particle model. 

We briefly describe the method together with a new regularization scheme in 
phase-space, which improves recovering the correct orbit structure in the models. We 
also mention some practical applications showing the power of the technique in inves- 
tigating the dynamics of galaxies. 



1. Introduction 

In the last decades, an increasing amount of high quality photometric and kinematic 
data for galaxies have become available. The dynamical state of the observed galaxy, 
however, cannot be directly inferred from observations, due to projection effects, and 
modelling is essential to learn about the distribution of stellar orbits and the total (i.e. 
due to luminous and dark matter) gravitational potential in the observed galaxy. There- 
fore, several methods to model the observational data and create "made-to-measure" 
models have been devised. 

For instance, assuming that all the integrals o f motion are known, one can fit obser- 
vations with parametrized distribution functions (| Deiongh e1|l984lll986l:[Bishorj|ll987l: 
Gerh ar3l9^lHunter & de Zeeuwfi99llCarollo et alJll995l:lKuiikedll995l :lMagorrian 
19951: iMerrittf 1 1 996k iDehnen & Gerhardlll993l:lMatthias & Gerhardlll999h. or solve the 



Jeans equations subject to the observational constraint s (Binney & Mamon 1982; Bin ney et al 
19901 : iMagorrian & Binnevil 19941 : lLokasll2002l : ICappeUarill2008h . " 

Another technique which is widely used is the Schwarzschild orbit superposi- 
tion method (ISchwarzschildl fl979[ 1 19931) : a large library of orbits is computed in a 
fixed potential, and then the weights of individual orbits are a djusted until the model 
matches the photometry and kinematics of the tar g et galaxy dRichstone & Tremaine 
1985 1: ICappellari et all 120021: iGebhardt et all 120031: IValluri et all 120041 : iThomas et al. 
20051: Ivan den Bosch & de Zeeuwll2010h . 

ISyer & Tremaind (Qj996) proposed an alternative, particle-based method. A mod- 
ified version sui t able f or modelling observational data with errors was designed by 
de Lorenzi et al.l (120071) and implemented in the pa rallel code NM AGIC . More recent 
imple mentations of this method can be found in iDehnenl (l2009h and lLong & Mad 
(120101) . The basic idea is to evolve a system of particles while slowly correcting the 
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individual weights of particles until the Af-body system reproduces the observational 
data. The method is very powerful, since no orbit library needs to be computed or 
stored, no symmetry restrictions need to be imposed, and the potential can be evolved 
self-consistently from the particles. Moreover, the algorithm properly accounts for ob- 
servational errors, and a great variety of observational data can be used simultaneously 
in the weight adaptation, including photometry, long-slit spectroscopic data, integral 
field absorption-line kinematics, and PNe velocities. 

So far, the particle made-to-measure m ethod has been used to investigate the dy- 
namics of the Milky Way's bulge and disk dBissantz et all 12004 ). an d the dark matter 
fraction and orbital structur e in the outer halos of elliptical galaxies ( de Lorenzi et all 
2008l . l2009l : lDas et aUl201ll) . 



2. Made-to-measure particle models of observational data 

An A?-body system of particles is trained to reproduce the observables of a target galaxy 
by maximizing the merit function 

F = -^x 2 +t*S (1) 



with r espect to the particle weights wt (see e.g. lSyer & Tremainelll996l : lde Lorenzi et all 



120071) . Such maximization strikes a compromise between the goodness of the fit (x 1 ) 
in terms of deviations between target and particle model observables, and a pseudo- 
entropy 

S = -2_jWilog(Wi/Wi), (2) 

i 

which serves the purpose of regularization. Since typically the number of particles is 
much higher than the number of noisy and sparse constraints on the particle model, the 
problem is intrinsically ill-conditioned, and some regularization is necessary. This is 
achieved by pushing particle weights towards a smooth distribution of predetermined 
priors w,-. 

The maximization of equation £[]) translates to a prescription for adapting the 
weights of the particles while they are evolved in the gravitational potential, which 
can be either fixed and known a priori, or time-varying and self-consistently computed. 

2.1. Regularizing particle models 

Traditionally, priors are set to vi>,- = wo = l/N (the "flat" priors in Bayesian statistics), 
and are kept constant during the modelling. Likewise, the individual weights of the 
initial particle model are usua lly set to 1 /N. 



However, as discussed in lMorganti & Gerhard! (l2012h . smoothing the weights glob- 



ally towards a set of preassigned, constant priors makes it difficult to construct made-to- 
measure galaxy models with both smooth and anisotropic DFs unless the initial particle 
model is already close to the target galaxy, i.e. the dynamics of the galaxy is already 
known. In this paper, we described and implemented in NMAGIC a new regularization 
method based on moving (i.e. not constant) priors which follow the smooth phase-space 
structures traced by the weight distribution, as the latter evolves to match the observa- 
tional data. The computation of such new individual priors is particularly easy when the 
integrals of motion are known, or at least one can compute approximately conserved 
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quantities from the particle orbits. In IMorganti & Gerhard (120121) we consider spheri- 
cal and axisymmetric cases, and show that the procedure facilitates recovering both a 
smoother and more accurate DF without erasing global phase-space gradients. 



3. Convergence to a theoretically unique solution 



It is a natural question whether the made-to-measure modelling technique can recover 
the target pr operties from observations i ndependently of the initial particle model. As 
we show in IMorganti & Gerhard! (120121). if a unique inversion of data to recover the 
underlying target DF exists dDejonghe & Merrittlll992r) . then it can actually be found 
via NMAGIC modelling. In Fig. Q] the quality to which the orbital anisotropy of the 
target galaxy can be recovered in such a case for different initial particle models can be 
appreciated. 
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Figure 1 . Recovery of the anisotropy parameter of a spherical target galaxy (full 
black line) from idealized line-of-sight velocity data. Radius is normalized by the 
model's half-light radius R e . The black lines, nearly on top of each other, show the 
final NMAGIC models obtained with our new regularization scheme using different 
initial particle models (corresponding grey lines). 



By contrast, lack or poor quality of data intro duce degeneracies in the dynamical 
modelling results; see IMorganti & Gerhard! (120121) for a quantitative analysis. 



4. Modelling the halos of elliptical galaxies 



The outer halos of elliptical galaxies are particularly interesting because they are dark 
matter-dominated, and because the orbital structure of stars there more strongly pre- 
serves the imprint of formation mechanisms, due to longer dynamical time scales. 

Dynamical models of these galaxies were constructed with NMAGIC, fitting a 
wide variety of photometric and kinematic data in a sequence of gravitational poten- 
tials for the dark matter halos. When modelling the interme diate luminosity ellipti- 
cals NGC 4697 and NGC 3379, Ide Lorenzi etaiI(l2008L [2009T) found that a variety of 
potentials ar e consistent with the da ta, showing that the well-known mass-anisotropy 
degeneracy dBinney & Mamon1ll982t) remains substantial when the velocity dispersion 
strongly decreases with radius, even with data out to more than 5R e . By contrast, the 
mass distrib ution in luminou s ellipticals with extended flat dispersion profiles is well- 
constrained ( Das et ai1l201ll) . 
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5. Conclusions 

The parallel code NMAGIC is a powerful tool to build made-to-measure galaxy models, 
which reproduce a wide variety of observational data. NMAGIC works by adapting the 
weights of an Af-body particle system until the target observables are well matched, 
subject to additional regularization constraints. 

We have briefly described a new, improved regularization scheme in phase-space, 
and shown how well the intrinsic properties of the target galaxy can be recovered inde- 
pendently of initial conditions, and how NMAGIC models can help us in learning about 
the orbital structure and gravitational potentials in galaxies. 
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